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Pushing NRQCD to the hmit 

Chris Stewart and Roman Koniuk ^ 
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Lattice NRQCD has proven successful in describing the physics of the upsilon system and B mesons, though 
some concerns arise when it is used in simulations of charm quarks. It is certainly possible that the NRQCD 
expansion is not converging fast enough at this scale. We present some preliminary results on the low-mass 
breakdown of NRQCD, in particular the behaviour of QQ and Qq spectra as the bare Q mass is decreased well 
below 1, with the aim of understanding more about the manifestation of this breakdown. 



1. Introduction 

Lattice NRQCD studies of heavy-quark sys- 
tems have been, on the whole, very successful. 
The predicted spectra for the Upsilon, J/^, B 
and D systems, amongst others, agree with the 
overall structure of the experimental spectra, and 
for the heavier systems, the agreement at finer 
scales, such as hyperfine splittings, is at the per- 
cent level 0]. 

Trottier's work on the charmonium spectrum 
m indicates that the 0(z;'^)-improved NRQCD 
action decreases the hyperfine splitting signifi- 
cantly over the O(w^) result, rather than taking it 
towards the experimental value. Our own results 
using similar parameters to Trottier confirm this 
result. This is an indication that NRQCD is not 
converging well for the charm quark. 

Yet it is interesting to note that many suc- 
cessful quark model predictions of the light-quark 
spectrum used a non-relativistic approximation 
for the light-quark dynamics B]. Here we have 
a highly relativistic system reasonably well de- 
scribed with a non-relativistic theory. 

Recently, Liu et al. |^ published work on a 
model QCD theory known as Valence QCD. In 
VQCD, all z-graphs are removed, and the authors 
were able to make some links between NR quark 
models and the role of z-graphs in standard QCD. 
We would like to examine the behaviour of the 
low-mass limit of NRQCD, to examine the nature 
of the inevitable breakdown of the NR expansion, 
and see if this can also be linked to the remarkable 



success of NR quark models. In this report we 
present results for the QQ and Qq spectra, as a 
function of decreasing heavy quark mass. 

2. Simulation details 

We used an C'(wS)-improved NRQCD Hamilto- 
nian, 
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where a A denotes the use of improved derivative 
operators, and E, B are components of the im- 
proved field tensor. The heavy-quark propagator 
was calculated using the evolution quation 
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For Qq mesons, we used the standard tadpole- 
improved SW action for the light quarks, with 
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Figure 1. Quarkonium ^Sq, "^Si and ^Pi states as 
a function of Mq . 



Figure 2. Quarkonium hyperfine and P — S split- 
tings as a function of Mq. 



uq = {Plaq)^'^ and Csiu = Uq . The light quark 
K — 0.143, roughly comparable to the strange 
quark mass. 

Gauge field configurations were created using a 
standard tadpole- and rectangle-improved action 
B. We used an ensemble of *** configurations 
of size 8^ X 24 at /? = 6.8, which gives a tadpole 
factor of uo = 0.854. 

3. Quarkonium 

Results for the ^S'o, '^5'i and ^Pi states of 
quarkonium for various values of Mq are shown 
in Figure 1. The stabelisation parameter n in 
Equation (1) was given values ranging from 2 to 
10 depending on the value of Mq. Meson opera- 
tors were smeared at source and sink, and their 
specific forms are detailed in P]. 

The ^5*0 states in Figure 1 are found using the 
kinetic definition of the meson mass, 



E-p — Eq 
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where Ep and Eq are the simulation energies for 



a meson with finite and zero momentum respec- 
tively. This method leads to the large error bars 
shown on the data, however the general trends 
are significant here, not the exact energies. 

The heaviest Mq corresponds to roughly half 
the bottom quark mass, while the charmonium 
ground state energy of ~ 3GeV corresponds to 
aAIo ~ 1.5 to 2. The lightest Mq gives a kinetic 
mass close to the energy of an ss meson, and so 
we can assume the quarks are well within the rel- 
ativistic regime at this stage. 

Note that the S and P states decrease quite 
linearly with Mq, until the bare mass drops sig- 
nificantly below 1, when they begin to drop quite 
suddenly. The hyperfine splitting increases sig- 
nificantly at low Mq, as might be expected from 
light-meson spectroscopy, but the S — P splitting 
decreases, against observation. These splittings, 
with the '^S'l state taken as the energy zero, are 
shown in Figure 2. 

These results for quarkonium indicate the 
NRQCD action is indeed beginning to have dif- 
ficulties at these low bare masses. The situation 
is worse for the heavy-light spectrum. Figure 3 
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Figure 3. Hyperfine splitting of Qq system as a 
function of A/q ■ 



shows the hyperfine spHttings for the Qq system, 
and note that for the lightest value of Mq, the 
splitting is negative — the '^Si and ^5*0 energies 
are in the wrong order. This is obviously a se- 
rious pathology. 

4. Discussion 

The results presented here are perhaps not 
surprising, given that we expect the Mq - 
dependence of the coefficients in the NRQCD 
action will ultimately lead to a breakdown as 
Mo decreases. However, we are interested in 
the exact nature of this breakdown, not to see 
if NRQCD is a viable alternative for simulat- 
ing light quarks, but to illuminate any connec- 
tions between NRQCD and the successes of non- 
relativistic quark models. 

Lewis and Woloshyn Q have shown that un- 
physically large contributions from certain terms 
in the NRQCD action could be removed by sub- 
tracting their expectation values from the Hamil- 
tonian. Their analysis consisted of a thorough 
systematic examination of the contributions of 



each correction term in the Hamiltonian. We sus- 
pect that the same style of analysis, applied to the 
low-mass breakdown of NRQCD, would fiush out 
the terms that suffer the worst pathologies, and 
may even suggest ways they may be strengthened. 

One important issue to be aware of, however, is 
that the coefficients in Equations (1) and (2) are 
usually given their tree-level value of 1 in NRQCD 
simulations. For heavy quarks, tadpole improve- 
ment of the gauge field links is usually sufficient 
to account for the most serious corrections be- 
yond tree level. The coefficients q in Equation 
(2) will, in general, have non-trivial as(Afo) de- 
pendence, and this will become more important 
as the quark mass is decreased. 

We would like to thank H. Trottier, R. Lewis, 
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tional Science and Engineering Research Council 
of Canada. 
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